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Abstract: We initiate the study of how to extend the correspondence between dimer 
models and (0+l)-dimensional cluster integrable systems to (1 + 1) and (2+l)-dimen- 
sional continuous integrable field theories, addressing various points that are necessary 
for achieving this goal. We first study how to glue and split two integrable systems, from 
the perspectives of the spectral curve, the resolution of the associated toric Calabi-Yau 
3-folds and Higgsing in quiver theories on D3-brane probes. We identify a continuous 
parameter controlling the decoupling between the components and present two comple- 
mentary methods for determining the dependence on this parameter of the dynamical 
variables of the integrable system. Interested in constructing systems with an infinite 
number of degrees of freedom, we study the combinatorics of integrable systems built 
up from a large number of elementary components, and introduce a toy model captur- 
ing important features expected to be present in a continuous reformulation of cluster 
integrable systems. 
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1. Introduction 

It has been more than a decade since the emergence of the technology of singularity 
resolution by D-branes, in the construction of supersymmetric gauge theories on their 
worldvolume. This connection between gauge theory dynamics and Calabi-Yau geom- 
etry has been a triumph of the AdS/CFT correspondence [|l], |], [|. Much progress has 
undoubtedly been made, especially for the largest class known to the correspondence, 



constituting almost all the known examples. This is the class of toric singularities 

11,1,0- 

For four-dimensional gauge theories arising from toric Calabi-Yau threefolds, the 
geometry is controlled by a planar toric diagram consisting of a convex, lattice polygon 
and its interior lattice points. On the other hand, it has by now become clear that the 
best way to understand the physics is through dimer models or, equivalently, periodic 



brane tilings on the plane || |], [10], p| . The complete encoding of the worldvolume 
quiver gauge theory in terms of these bi-partite graphs drawn on a torus has led to 
both the simplification of existing problems as well as development of new ideas in 
a plethora of directions ranging from physics to mathematics. The computation of 
moduli spaces has been trivialized, quantities such as R-charges or procedures such 



as Higgsing or Seiberg duality |12| have found new graphical realizations. In parallel, 
geometric resolution by combinatorics of perfect matchings or number-theoretic issues 
such as the field extension of the defining torus have also provided a fruitful dialogue. 

Last year, the path to a new direction was illuminated in the mathematics literature 
by [[H|, where a correspondence between dimer models and certain integrable models 
was found. The correspondence associates any dimer model on a torus to a (0 + 1)- 
dimensional quantum integrable system - which they dub cluster integrable system - 
whose phase space contains the moduli space of line bundles with connections on T and 
whose Hamiltonian and Casimir operators are given by the partition function of the 
dimer model. 

This connection was readily exploited in our current context of toric gauge theories 
in |TJ], I5J, |IE |. Explicit integrable systems can be thus established for the myriad of 
theories in the catalogue (cf. |C7L 0) of brane tilings. 

Fortified by the strengthening of this correspondence in various guises, we propose 
to extend it in yet another direction, investigating whether a certain continuous limit 
exists. Indeed, our dimer models are finite graphs, thereby giving a finite number of 
faces in the fundamental domain of the torus. Loops around these faces, together with 
loops around the two fundamental directions of the torus, provide a parametrization of 
the Poisson manifold of the integrable system. Nevertheless, one could envision a toric 
diagram with an infinite number of lattice points corresponding to some singularity of 
infinite order - obtained for example by taking the limit of an Abelian orbifold with 
the quotient group taken to infinite order. In this situation, we shall be dealing with 
an integrable system with an infinite number of degrees of freedom. Can we construct 
in this way continuous integrable field theories in (1+1) or even (2+1) dimensions? If 
so, can dimer models provide new tools or useful perspectives? 

In order to initiate this study we need to understand how to "glue" and "split" 
two integrable systems in our present framework. The procedure of splitting or gluing 



toric diagrams, in relation to geometric resolutions, is well-known. An efficient method 
for studying general resolutions exploiting dimer models was introduced fll9| . In this 



paper we will first systematically study how this splitting/gluing is mapped to the 
integrable-system variables, and will demonstrate how this manifests in the spectral 
curve (defined as the zero locus of the Newton polynomial) associated to the toric 
diagram. We will also discuss the combinatorics of contributions to conserved charges 
in the limit in which a large number of building blocks are pasted. Finally, we will 
introduce a toy model that exhibits some of the desired features of a reformulation of 
cluster integrable systems in the continuous limit. 

The organization of the paper is as follows. In Section |2] we review the correspon- 
dence between dimer models and integrable systems. Section ^| discusses the gluing 
and splitting of spectral surfaces and its relation to the resolution of toric singulari- 
ties. In |] we identify a continuous parameter controlling the separation of individual 
components. We relate this parameter to expectation values in the underlying quiver 
gauge theory and introduce two methods for finding the dependence of the integrable 
system on it. These ideas are illustrated with explicit examples in Section |5|. In Section 
D, we investigate the combination of a large number of building blocks and introduce 
a simple toy model that exhibits some properties we expect in the continuous limit of 
cluster integrable systems. We conclude in Section [7[ 



2. Dimer Models and Cluster Integrable Systems 

A remarkable correspondence linking dimer models to an infinite class of integrable 



systems, denoted cluster integrable systems, was recently introduced in p3| . We now 
provide a brief review of the correspondence. 

The Poisson manifold of the integrable system is parametrized by oriented loops 
on the brane tiling. Cycles going clockwise around each face Wi (i = 1, . . . , N g , with 
N g the number of gauge groups in the quiver) and the cycles z\ and z 2 wrapping the 
two directions of the 2-torus provide one possible basis for loops. The ^-variables 
are subject to the constraint rL=i w i = 1> which sometimes can be used to simplify 
expressions. 

The Poisson brackets are given by 

{wi,Wj} = e WhWj WiWj 

{ Zl ,z 2 } =l + e 2li22 (2.1) 



where e XyV is the number of edges on which the x and y loops overlap, with orientation. 
Then, e Wi)W . is simply the antisymmetric oriented adjacency matrix that counts the 
number of arrows between gauge groups in the quiver dual to the brane tiling. 

The integrable system can be quantized replacing the Poisson brackets by a q- 
deformed algebra, which takes the form 

XiXj = q^XjXi , (2.2) 

where Xi = e Xi , q = e~ t2nh and n^ = {xj, Xj}/(xiXj). 

Every perfect matching is associated with a point in a toric diagram [BJ. When 
more than one perfect matching corresponds to the same point, we must add all their 
contributions. As we have just reviewed, the natural variables of cluster integrable 
systems are loops. We can translate any perfect matching into a closed loop on the 
tiling by subtracting a reference perfect matching. 

In |JT5| , it was shown that the commutators defined by ( |2.2| ) and ( |2.1| ) give rise to 
(0 + l)-dimensional quantum integrable system, whose conserved charges are: 

• Casimirs: they commute with everything and are given by the ratio between 
contributions associated to consecutive points on the boundary of the toric dia- 
gram. 

• Hamiltonians: they commute with each other and correspond to the internal 
points in the toric diagram. 

The toric diagram of the Calabi-Yau 3-fold associated to the dimer model gives rise 
to a Riemann surface of genus equaling to the number of internal points. The latter is 
given by the zero locus of the Newton or characteristic polynomial 

D = {{v l7 v 2 )i} ~> P(z h z 2 ) = Y^Cizfz? , (2.3) 

i 

where D is the planar toric diagram, specified by lattice points (v 1,1*2)1 and zi,z 2 € <D 
are complex coordinates. The coefficients q are functions of the w-variables. This 
Riemann surface is indeed the spectral curve of the integrable system. 

The full Poisson manifold of the integrable system is obtained by gluing differ- 
ent patches via cluster transformations, equivalently Seiberg duality in the associated 
quiver gauge theories. 



3. Gluing and Splitting 

In this section we discuss the decomposition of a spectral curve into pieces and the 
reverse procedure of gluing spectral curves. We also consider another context in which 
this process arises, the desingularization, or Higgsing, of a Calabi-Yau space. The 
intimate connection between the decomposition of integrable systems and Higgsing 
will be the topic of forthcoming sections. 

3.1 Spectral Curves 

Let us consider the splitting process 



Si 



(3.1) 



where £ is the "parent" spectral curve and S 1)2 are the two daughters. In this process we 
elongate certain throats of the E until it breaks into two pieces. Since the spectral curve 
is a thickening of the (p, g)-web p(J ^TJ, which is the graph dual to the toric diagram, 



the connecting throats are dual to segments joining points of the toric diagram along 
the boundary between the daughters, as shown in Figure |l|. 





Figure 1: Splitting of a Riemann surface into two daughters. 



The separation between components of the spectral curve is achieved by tuning 
the coefficients q in the characteristic polynomial fl2.3|) . In the integrable system, these 
coefficients are functions of the w;-variables. In the limit of large distance between 
components, the q's scale differently with respect to the separation and develop a 
hierarchical structure. 



Given a decomposition of the spectral curve into two pieces, it is natural to expect 
that in the limit of infinite separation the original integrable system reduces to the sum 
of the integrable systems associated to the two components. In Sections |] and ||, we 
provide a detailed explanation of how this intuition is realized. 

Taking the splitting process to an extreme we obtain a decomposition of any Rie- 
mann surface into a collection of trinions, i.e. spheres with three punctures. Any such 
decomposition is in one-to-one correspondence with triangulations of the toric diagram. 
Each triangle gives rise to a trinion. Their number is thus equal to twice the area of 
the toric diagram, which in turn is equal to the number of gauge groups in the asso- 
ciated quiver gauge theory. These decompositions allow us to see the full integrable 
system continuously emerge from the combination of trivial integrable systems asso- 
ciated to trinions. Figure |2] shows a possible triangulation of a toric diagram and its 
corresponding trinion decomposition. 




Vr 



Figure 2: A general triangulation of the toric diagram and the corresponding trinion de- 
composition of the spectral curve. 

A standard method in real ("tropical") geometry to visualize Riemann surfaces, 
which will prove useful in later sections, is the so-called amoeba projection: 

A : (z 1 ,z 2 ) •->• (log | zi |, log|,z 2 |). (3.2) 

The amoeba can be thought of as a thickening of the graph-dual of the toric diagram; i.e., 
we can draw the (p, q)-web from D and this will constitute the "spine" of the amoeba. 
In the actual plot, the "tentacles" which tend to infinity will have their directions given 
by the (p, g)-vectors which are normal to the toric diagram. 

As explained in the introduction, what we ultimately wish to study is the integrable 
system that emerges in the continuous limit where we glue a countably infinite number 



of toric sub- diagrams, or equivalently, spectral curves. We conjecture that, depending 
on how we assemble these building blocks, the integrable systems with an infinite 
number of degrees of freedom that are generated by this procedure are (1 + 1) or 
(2 + l)-dimensional integrable field theories. Figure || shows a number of elementary 
spectral curves glued to generate a (1 + l)-dimensional theory. 

Figure 3: Combining an infinite number of elementary spectral curves to generate a (1 + 1)- 
dimensional theory. By performing a similar gluing along the vertical direction we expect to 
generate a (2 + l)-dimensional theory. 



3.2 Partial Resolution of Calabi-Yau Singularities and Higgsing 

The splitting and gluing of Riemann surfaces of the type discussed in Section [D] also 
play an important role in the context of partial resolutions of toric singular Calabi- 
Yau 3-folds. We now review this process from geometric, gauge theoretic and dimer 
model perspectives. In Section [| we elaborate on the intimate connections with the 
(de) composition of integrable systems. 

Geometrically, the partial resolution of a toric singularity corresponds to the process 
illustrated in Figure [L| one takes the toric diagram and divides it into components 
1 . The resulting components must correspond to well-behaved toric diagrams, which 
constrain them to be convex. For concreteness, we will focus on the case in which 
we split the toric diagram into two parts, to which we shall affectionately refer as the 
parent with two daughters. It is possible to deal with more components by iteration of 
this procedure. 2 

Let us now describe the resolution from the perspective of the quiver theory on the 
worldvolume of D3-brane probes. This is of course standard technology dating back 
to the early days of studying quiver gauge theories from toric Calabi-Yau singularities 
IS H H> 0- The starting point is a set of N = n\ + n 2 D3-branes on the parent 
singularity. All chiral fields in the parent quiver are (ni + n 2 ) x (ni + n 2 ) matrices. The 



1 In the dual cone picture of the toric variety, this is the process of stellar division p2[ . 
2 Not all decompositions can be reduced to a sequence of binary splittings. A necessary condition 
is that, at each step, the toric diagrams of the daughters are convex. 



parent singularity is then resolved into two daughter singularities containing n\ and n 2 
D3-branes, respectively. As a result, we obtain two decoupled quiver gauge theories, 
whose gauge group ranks are given by n\ and n 2 - From a gauge theory viewpoint, 
this resolution corresponds to turning on non-zero vacuum expectation values (vevs) 
for some block sub-matrices in the scalar components of these fields. Fields charged 
under gauge groups in both quivers have masses controlled by the expectation values 
and decouple from the low energy theory. 

A pictorial representation of this process is given in Figure [|| One could envisage, 
of course, the reverse process of gluing to produce a more singular parent. This should 
correspond to an un-Higgsing mechanism (q. v. f23|). We will illustrate these ideas 
with explicit examples in Section [|. 




N=n l4 .n 



l+"2 



Figure 4: The resolution of the parent singularity with parallel coincident N = ni+ri2 D3- 
branes (on the left) results in two daughter theories with n\ and n 2 D3-brane respectively. 
The blue cone signifies the Calabi-Yau singularity at the tip of which the D3-branes sit. 



For our purposes, it is sufficient to focus on the simple case in which rti = n 2 
We focus on diagonal vevs of the form 



w 



x 



(i) 







X 



(2) 
ij 



(3.3) 



for the field X^ in the parent theory. We will stick to this case throughout the paper. 
Resolutions generically involve turning on several non-zero expectation values simulta- 
neously. We restrict to the case in which all non-zero vevs have the same magnitude, 
which will turn out to control the distance between the daughter singularities. 3 The 
acquisition of vevs in this fashion will split the parent theory into its two daughters. 



3 Theories with different vevs give rise to multiple energy scales. If these scales are hierarchically 
separated, the Higgsing process can be studied sequentially. 



3.2.1 Resolution in the Dimer Model 

One of the greatest computational challenges to the resolution of toric singularities by 
D3-branes was the identification of which fields in the parent acquire non-zero vevs 
H 0. This issue was resolved by the dimer model representation of toric quiver gauge 



theories || |9], |TT| . Dimer models are extremely useful for identifying the non-zero vevs 
that are necessary in order to achieve a given resolution. An elegant description of 
general partial resolutions, exploiting the map between the dimer model and a tiling 



of the spectral curve, was introduced in |1|. We now briefly review this procedure. 

Like perfect matchings, zig-zag paths play a prominent role in connecting dimer 
models to geometry. They are defined as paths that alternate between turning max- 
imally right and maximally left at consecutive nodes in the bane tiling. Each edge, 
then, has exactly two oppositely oriented zig-zag paths, criss-crossing before heading 
to nodes of opposite color, weaving an intertwined pattern on the torus. For consistent 
gauge theories, these zig-zag paths never intersect themselves and form closed loops 
wrapping (p, g)-cycles on the torus. 

The untwisting map is an operation on zig-zag paths that exchanges the criss-cross 
and turns a brane tiling, which by construction lives on a 2-torus, into a tiling of 
the spectral curve E (which we recall is the zero locus of the Newton polynomial for 
the given toric diagram) and vice versa. We refer the reader to |T]J for a detailed 
explanation of the untwisting map whose effect on zig-zag paths of both T 2 and E is 
summarized below. 

T 2 E 



zig-zag path o face = puncture 

face = gauge group 4-> zig-zag path 

Starting from the parent spectral curve E, we elongate one or several internal 
throats that connect the two daughters, Ei and E 2 . The daughters then decouple 
in the limit in which these throats become infinitely long. From the viewpoint of the 
daughters, these throats become new external legs, i.e., new punctures. The appearance 
of new punctures can easily be implemented in terms of brane tilings. We consider one 
copy of the original brane tiling on T 2 for each of the two daughters. On each copy, 
we draw the zig-zag paths associated to the original punctures that will end up on the 
corresponding component. Next, we introduce the paths which are the complement 
to these zig-zag paths in the original set. These new paths correspond to the new 
punctures that are generated in the splitting process. 

In order for the new paths to become actual zig-zag paths, some edges must be 
removed from the daughter tilings. The bifundamentals on the tiling that do not have 



any paths running over them are removed. These are precisely the ones that acquire 
non-zero vevs in the Higgsing. Generically, different edges are removed from the two 
tilings associated to the daughters. This is the manifestation, in dimer language, of 
the matrix vevs in ( |3.3| ). In Section |^, we present explicit examples illustrating this 
procedure. Bivalent nodes might be generated when removing edges. They correspond 
to massive fields that can be integrated out [BJ . From the perspective of the daughter 
integrable systems, the generation of new punctures corresponds to the appearance of 
new Casimir operators. 

When Higgsing, the rule for removing perfect matchings is simple: every perfect 
matching containing an edge corresponding to a field with a non-zero vev, must be 
eliminated. The mapping of the dimer model into a tiling of the spectral curve pro- 
vides an interesting alternative perspective on the disappearance of perfect matchings. 
Differences of perfect matchings are translated into 1-cycles with appropriate homology 
on E. The difference between two adjacent perfect matchings, say p and p', in the toric 
diagram simply corresponds to a 1-cycle wrapped around the associated elongation of 
E. Once E has been split, certain 1-cycles will no longer be able to exist. If the differ- 
ence between two perfect matchings is contained in any of the sub-dimers, both perfect 
matchings will survive in the corresponding component of the daughter singularity. On 
the other hand, if the difference is not contained in any sub-dimer, one or both of the 
perfect matchings will not survive the Higgsing process. 

In general, the set of non-zero expectation values resulting in a given decomposition, 
or, equivalently, the set of removed edges from the daughter brane tilings, is not unique. 
In Section [5], we will discuss the issue of multiple solutions for an explicit example in 
detail. 



4. A Continuous Control Parameter 

As we have already mentioned, the splitting of the spectral curve follows from certain 
hierarchies between the coefficients in the characteristic polynomial. These hierarchies 
are controlled by a continuous parameter, which we will denote A; in quiver language 
it is connected to the non-zero expectation values of bifundamental fields. 

In this section, we introduce two complementary approaches for determining the 
precise dependence of A on the coefficients of P(zi, z 2 ) that achieves a given E — > EX+E2 
decomposition. 



10 



4.1 Scalings from VEVs 



(<0\ 



Without loss of generality, we can identify A with the non-zero vevs, i.e. (JQ- ) = A. 



The coefficients in P(z\, z 2 ) are polynomials in the Wi variables, corresponding to closed 



loops with vanishing homology on the brane tiling. Following |24) (see also |25] for 



applications of this idea), we can write any loop in terms of edges in the brane tiling as 

J* a(w 4+ i, bj) 

where the product runs over the contour 7 and bj and Wj denote black and white nodes. 
Here, X(wj, bj) and X(wj +1 , bj) are bifundamental fields, in which we explicitly indicate 
the tiling nodes connected by the corresponding edge when going around 7 instead of 
employing the usual notation with subindices for the gauge groups under which they 
are charged. Going back and forth between the two notations is straightforward. Fur- 
thermore, we remind the reader that the variables defined Section [^ for our integrable 
system, in terms of the notation in ( |4.1| ), are 

Wj = v(j w .) ; Zi = v(~t z .) j = l,...,N g , i = 1, 2 . (4.2) 

We can thus obtain the A-scaling of any loop by plugging the expectation values 
into (|4. 1|) . The conclusion is that we obtain the following factors 4 

Scaling A for each (AT(wj, bj)) = A 
Scaling A" l for each (X(w i+1 , b*)) = A 

Applying fl4.3|) directly to the Wi cycles, we can rephrase it in terms of vevs for 
fields transforming in the fundamental or antifundamental representation of the corre- 
sponding gauge group. We obtain 

within Wi cycle : A for each (Xji) = A , . 

A" 1 for each (xj) = A '' 

where we have changed the notation and subindices indicate gauge groups connected 
by bifundamentals. We show some examples in Figure |5|. 

Notice that we have not included color indices in (|4.1| ). Similarly, we have not 
specified on which of the two daughter components, which have different sets of non- 



zero vevs, the scalings in (|4.3|) and (|4.4|) must be calculated. However, on physical 



4 Exchanging black and white nodes, an operation that has no effect on the physics, exchanges A 
with A -1 . This operation does not affect the scalings of individual contributions to conserved charges 
though, since it also inverts the orientation of all paths. 
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A A- 1 A 

Figure 5: Examples of A scalings of a Wi cycle, which we have drawn as the cyclic arrow in 
blue. We indicate the edge associated with a field with a non-zero vev in red. 

grounds, it is clear that determining the scaling in any of the two daughters leads to 
the same result. We shall present explicit examples in Section [5]. 

We have explained how to determine the A of any loop in the tiling arising from a 
general decomposition of the Riemann surface. It is important to emphasize that the 
A-independent part of these loops can be completely arbitrary in size, but in a way 
that is uncorrelated with the splitting of S. 

4.2 Scalings from the Spectral Curve 

Let us now introduce an alternative method for determining the A scalings directly 
associated with a decomposition of the spectral curve in an algorithmic fashion, wherein 
we obtain a set of explicit Diophantine inequalities which needs to be solved. Suppose 
the parent spectral curve, corresponding to the toric diagram Dp = {^1^2} (with i 
indexing the nodes) can be written as 

p M=EEn^ W ; (4 - 5) 

i \k=i j=i j 

here, we have explicitly written the coefficients as sums over monomials in the w- 
variables. Indeed, the number of terms, which we index by k, of monomials for the i-th 
node is the number of perfect matchings pi for that node. Finally, each monomial is a 
product of w- variables, indexed by j, raised to integers a\ -. 5 

Now, when we Higgs, the coefficients - as polynomials in ^-variables - separate 
into the sum of two types of terms: those which survive the Higgsing, which we shall 
separate into the first s, terms for the i-th node, and those which do not, which are the 
remaining p t — Sj terms: 

/ Si N G Pi N a \ 

p&> *>) = E E n < k " + E n 4" #4 ; (4-6) 

i \k=l j=l k= Si +lj=l J 



> by using ]J ,-= 



5 It is always possible if desired, although not necessary, to bring the monomials to a form in which 

4j 
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as an emphasis, we have relabeled the powers of the monomials which do not survive 
as ft d . 

Now, let us introduce the scalings (we can take, without loss of generality, the 
powers kj to be integers): 

Wj — ► A Kj , Kj e TL (4.7) 

and substitute back into ( |4.6|) to give 



Nr 



JV G G 

fc=l k=Si+l 

It is now clear what has to occur: for the terms which survive, they must have 
the same order in A and those which do not, must have strictly less order. In other 
words, for each i, we must have Ylj=i K j a \j = S?=i K j a 2,j = • • • = Si=i K j a li,j — C- 
This C must be strictly greater than each of Y^ij=i K jPk,j for fc = Sj + 1, . . . ,pj. Finally, 
we recall that JT-j^Wj = 1, so that Ylj=i K j = 0- ^he terms that do not survive are 
precisely those associated to perfect matchings that are removed by the Higgsing. We 
can use this identification to determine a priori which terms must have a relative A 
suppression. 

In summary, we have the following set of Diophantine inequalities in k: for each 
i — 1, 2, . . . , n where n is the number of nodes in the toric diagram, 

N G 

C = y^ K j a k,j f° r all fc = 1, . . . , Sj 

i=i 
N G 

c > y^ KjPij f° r eacn k = Si + 1, . . . ,pi 

i=i 
n g 



5. Explicit Examples 

Having abstractly discussed how the splitting should work by obtaining the weights of 
the variables either from the acquisition of vevs in the dimer or from the coefficients 
in the spectral curve, we can now illustrate our proposal in detail with some explicit 
examples. In this section, we will first study the cone over the double zeroth Hirzebruch 
surface and then the space Y 4,0 . These will give ample demonstration of our technique. 
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5.1 Double F 

Let us consider a Z 2 orbifold of F , whose toric diagram is shown in Figure |7|, to which 
we refer as the double Fq theory. Figure ^| shows the corresponding brane tiling. We see 



%> 



-*•£> 



\U- 



-5' 



K> 



<> 



Qtt*-^* 



K> 



-6 



-6 



Figure 6: Brane tiling for the double -Fq theory. 



that there are 8 the gauge groups, twice that of F , and 16 bifundamental fields which 
we shall denote as Xij in standard nomenclature, signifying the field corresponding 
to the edge bounding face i and face j in Figure []. The superpotential terms are all 
quartic and can also be instantly read off from the figure. We wish to consider the 
decomposition of this geometry into two copies of Fq, as shown in Figure [7[ 




Figure 7: Toric diagram for the double Fq. A red line indicated how we split it into two 
components. 



Using the ideas of [19], which were summarized in Section p. 2.1 , we conclude there 
are four possible sets of expectation values leading to the same desired decomposition 
of the geometry. They are: 



Higgsing 1: 



Higgsing 3: 



y(l) yW yM y« 
^Ml 5^12 5^85 5^56 



Y (2) Y (2) 
^83 ' ^32 



Y (2) Y (2) 

, yv 47 , w 76 



v (l) Y (i) Y (i) Y (l) 
^27 5^78 '^63 5^34 
(2) v (2) v (2) v (2) 



yyt) Y \z, ) y^j y- 

^25 ?^54 5^61 i-^ 



Higgsing 2: 



Higgsing 4: 



18 



y(l) vr(l) y(l) y(l) 
^34 5^-41 )^78 5^85 
y (2) y (2) y (2) Y {2) 
^32 5^25 5^76 5^61 



y(l) y« yW yW 
^-12 )^27 '^56 '^63 
y (2) y (2) y (2) y (2) 

^18 '^83 5^54 5^47 



(5.1) 
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We have separated, in the above, the 8 fields of the parent theory that get a non-zero vev 
into its two daughters, which following the notation in ( |3.3| ) we denote by superscripts 
(1) and (2) respectively. We see that each daughter contains four fields with non-zero 
vevs. From these non-zero vevs and the rule prescribed in ( [4.4] ), we determine the 
weights of our Wi variables for the four Higgsings: 



Higgsing 


Wx 


U>2 


w 3 


W{ 


w 5 


w 6 


wr 


w 8 


1 


1 


A 


1 


A- 1 


1 


A 


1 


A- 1 


2 


A 


1 


A" 1 


1 


A 


1 


A" 1 


1 


3 


1 


A" 1 


1 


A 


1 


A" 1 


1 


A 


4 


A- 1 


1 


A 


1 


A- 1 


1 


A 


1 



(5.2) 



What is happening in the field theory, as graphically depicted by the dimer, is 
shown in Figure [8]. In each of the four Higgsings, we separate the parent dimer model 
into the complementary dimers of the two daughters, (1) on the left and (2) on the 
right. In order to facilitate comparison with the original parent tiling, we have not 
integrated out massive fields. If we do so, we obtain the square lattice characteristic of 
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Higgsing 1 



Higgsing 2 
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Higgsing 3 



Higgsing 4 



Figure 8: Dimer models for the four Higgsings of the double Fq theory into two daughter 
Fq theories. 

The weights associated to different Higgsings are simply related by an overall shift, 
which follows from the fact that the Higgsed tilings are also connected by shifts and 
rotations as shown in Figure S. 
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The Integrable System 

We now investigate the effect of continuously splitting the spectral curve, equivalently 
Higgsing, on the integrable system. In this process some contributions to conserved 
charges, those associated perfect matchings removed by the Higgsing, are continuously 
suppressed until the theory reduces to two decoupled integrable systems. 

We can also determine the A-scaling of w-variables by analyzing the behavior of 
coefficients of the spectral curve. Indeed, we can find the weights in (|5.2|) independently 



by the procedure outlined in Section |4~2| . Reassuringly, the results of using this method 
agrees in all examples with the ones obtained from ( |4.4| ). 

The tables below present the information in the integrable system for the double 
Fq theory. The (r^i,^) row refers to the coefficient for the z^z^ 2 term in the Newton 
polynomial. For each of the Higgsings, we underline the contributions that survive in 
the A — > oo limit. The third column shows the surviving leading A dependence of the 
coefficients in the characteristic polynomial. 



Higgsing 1: Let us begin with Higgsing 1 and discuss it in a little more detail. For 
each node in the toric diagram, here given as an integer 2-vector, we can specify the 
coefficient in terms of the w-variables within the spectral curve using the technique of 
|15| , as reviewed in Section [| The number of monomials will correspond to the number 
of perfect matching for the node. Let us take the (1,1) point of the toric diagram as an 
example, which is an internal node with 8 perfect matchings; the term in the spectral 
curve will be 



U7lU7 3 U> 4 + U7 3 U7 4 + W3 + 1 + W Q 1 + W 5 X Wq 1 W % * + W1W2W3W4 + ^5 1 Wq l J 



Z1Z2 , 



where we have underlined the terms which survive the Higgsing. The ones that do not 

survive can be immediately determined, they are the ones containing edges correspond- 
ing to fields with non-zero vevs. 

In terms of ( |4.9| ), this means the weights Wj — > A k] must be such that k$ = = 
— k 5 — k e — ks = £4 + ^2 + ^3 + ^4 coming from the underlined terms and that they must 
all be strictly greater than any of {£4 + k 2 + k 3 , £; 3 + £; 4 ,—ke , —k 5 — k 6 } coming from 
the non-underlined terms. This thus constitutes one of the inequalities. We do this for 
each of the 8 points in the toric diagram and combine all these relations, supplementing 

8 

by the inequality that Yl kj = 0> an d solve the resulting system over the integers. We 



will find precisely the first row of the solution table in (|5.2|) . In the table below, we will 



16 



also include, for reference, the final leading order weight for the surviving terms. For 
the (1,1) term above, this is just 0, whence the entry A = 1 in the third column. These 
results are in full agreement with those derived using (| 



and (f4.4| ). In the examples 
that follow, we have independently determined the A-scalings using both methods and 
confirmed their agreement. 



(ni,n 2 ) 


Loops 




(0,0) 


1 


1 


(1,0) 


1 + WiW 2 W 5 W 6 


A 2 


(2,0) 


WiW 2 W 5 W 6 


A 2 


(1,1) 


WiW^Wa + W3W4 + W3 + 1 + Wq 1 
+W^ 1 Wq 1 W 8 1 + W1W2W3W4 + W^Wq 1 


1 


(2,1) 


WiW 3 + WiWq 1 + W$ + Wi + W^Wg 1 
+Wq 1 Wy 1 W% 1 + WiW 2 W 3 W 5 + WiW 2 W 3 


A 


(1,2) 


WsW^W^Wq 1 


A" 2 


(2,2) 


WiWzW^Ws 1 + WzW^W^Wq 1 


1 


(3,2) 


WiWsW^Wg 1 


1 



(5.3) 



Higgsing 2: We can now perform a similar analysis for the second Higgsing and obtain: 



(ni,n 2 ) 


Loops 




(0,0) 


1 


1 


(1,0) 


1 + WiW 2 W 5 W e 


A 2 


(2,0) 


WiW 2 W 5 W 6 


A 2 


(1,1) 


W1W3W4 + W3W4 + W 3 + 1 + Wq 1 
+W^ 1 Wq 1 W 8 1 + W1W2W3W4 + W^Wq 1 


1 


(2,1) 


wiw 3 + WiWg 1 + Wg 1 + wx + w 6 1 w s 1 

+Wq 1 Wy 1 W 8 ~ 1 + WiW 2 W 3 W 5 + WiW 2 W 3 


A 


(1,2) 


WzWaW^Wq 1 


A- 2 


(2,2) 


WiW-sWiWs 1 + WzW^W^Wz 1 


1 


(3,2) 


WiW 3 Wq 1 W & 1 


1 



(5.4) 
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Higgsing 3: So too we can now study the third Higgsing, confirming our results: 



(ni,n 2 ) 


Loops 




(0,0) 


1 


1 


(1,0) 


1 + W 1 W 2 W 5 W6 


1 


(2,0) 


WiW 2 W 5 W 6 


A" 2 


(1,1) 


WiW^Wi + W3W4, + W3 + I+ Wq 1 
+W^ 1 W 6 1 W S 1 + W 1 W 2 W 3 W 4 + W^Wq 1 


A 


(2,1) 


W1W3 + WlWg 1 + W$ + W\ + Wq Wg 


1 


+Wq 1 Wj 1 W^ 1 + W1W2W3W5 + W1W2W3 


(1,2) 


w^w^w^w^ 1 


A 2 


(2,2) 


WxWzW^Wq 1 + WsW^W^Ws 1 


A 2 


(3,2) 


VUxW^W^Ws 1 


1 



(5.5) 



Higgsing 4: Finally, we complete the story with the last Higgsing: 



(ni,n 2 ) 


Loops 




(0,0) 


1 


1 


(1,0) 


1 + WiW 2 W 5 Wq 


1 


(2,0) 


WiW 2 W 5 W 6 


A" 2 


(1,1) 


Wi%V 3 W A + W 3 W 4 + W 3 + 1 + Wfi 1 
+W^ 1 W 6 1 Wg 1 + W1W2W3W4 + W^Wq 1 


A 


(2,1) 


W1W3 + WiWg 1 + w$ l + Wi + W^Wg 1 


1 


+Wq 1 Wj 1 W^ 1 + WiW 2 W 3 W 5 + WiW 2 W 3 


(1,2) 


W 3 WiW^ l W G l 


A 2 


(2,2) 


WiWzWiWq 1 + WsW^W^W^ 1 


A 2 


(3,2) 


w 1 w 3 Wq 1 w 8 1 


1 



(5.6) 



Amoeba Projections 



As reviewed in Section ITT, amceba plots provide a simple way of visualizing the spectral 
curves. Thus, we can draw such projections of the spectral curve of the parent for some 
large value of A, as an additional check that the scalings we obtained indeed give rise 
to elongations leading to the desired splitting of the spectral curve. 

Let the A-weights for each node in the toric diagram be given by the third column 
in the tables above. The leading behavior in A coincides for Higgsings 1 and 2 and 



for Higgsings 3 and 4. Furthermore, the two pairs are related to each other by a 180° 
rotation, as shown in Figure pi. 





A^ 



1 and 2 



3 and 4 



Figure 9: The corresponding weights of the coefficients in the Newton polynomial for Higgs- 
ings 1 and 2 are connected to those of Higgsings 3 and 4 by a 180° rotation of the corresponding 
toric diagram. 



We see that, in perfect agreement, the amoeba projections exhibit the correspond- 
ing behavior. Figure [L0] shows the amcebas for the four Higgsings for large A. Indeed, 




Higgsings 1 and 2 




-12- 



Higgsings 3 and 4 



Figure 

e 5 . The 
in their 



10: Amoeba plots for the 4 possible Higgsings with A set to the numerical value of 
patchy appearance of these and subsequent amoeba plots, with some missing points 
interior, is due to the fact that we determine them numerically. 



the amcebas for the four models coincide, up to a trivial shift on the (x, y) plane. This 
results from the simple relation between their scalings as given by Figure §. Further- 
more, we see that the thin spine in the center controls precisely the spliting of the 
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double F into her two daughter F theories. The holes in the spectral curve associated 
to internal points in the toric diagram have zero size in Figure |10[ This is due to the 
particular choice of coefficients in the characteristic polynomial. These coefficients can 
be varied at will without modifying their A-scaling and hence preserving the splitting. 



1,0 



5.2 Y 

Fortified by the consistency our story for double F , let us move onto another non- 
trivial example. We now consider F 4,0 , whose corresponding integrable system has 



been worked out in |15[ . We remind the reader of the dimer model in Figure [II]. This 
is a theory with 8 gauge group factors, quartic superpotential terms and 16 fields, which 
we suggestively label as V^, V^ and H^ in accordance with their vertical and horizontal 



orientation. The toric diagram is shown in Figure 12, with 7 nodes. 
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Figure 11: Brane tiling for Y ' . 



We now study the two splittings depicted in Figure [T2]. We will refer to them as 
Higgsings 1 and 2. Specifically, the acquisition of vevs is as follows: 



Higgsing 1 : H&\H&\H$,H$ 

rr(2) „(2) rr(2) rr(2) 
-"34 ) -"67 ' ""12 ' -"85 



Higgsing 2 : V 5 



(i) 

51 



y{l) yW t>(1) 



(2) T ~H2) T/ (2) T/ (2) 
51 ' v 73 ' v 2& ' v 48 



(5.7) 



We remark that the black and white nodes are exchanged with respect to [15 |. 
This is of course only a matter of convention. As in the previous example, the choice of 
vevs leading to each splitting is not unique. Having already illustrated this possibility 
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Higgsing 1 



Higgsing 2 



Figure 12: Toric diagram for Y 4 ' , showing the two splittings that we will investigate. The 
splitting is indicated by the red line. 



with the splitting of double F$ to two i*o's, we focus on the vevs given in ( |5.7|) for 
convenience. 

Once again, we determine the A-scalings associated to both decompositions using 
the methods of Section [|, i.e. both the rule in Q4.4|) and the algorithm in Section [O 
which yield identical results. We obtain: 



(5.* 



Higgsing 1: We follow the notation in the previous examples, write the coefficients 
in terms of the w-monomials and underline the terms which survive. We tabulate this 
for each node, and in the third column write the overall leading behavior in A for the 
terms which survive. 



Higgsing 


W\ 


w 2 


w 3 


tt>4 


w 5 


Wq 


W'j 


w s 


b 


A- 1 


A 


A- 1 


A 


A 


A- 1 


A 


A- 1 
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A 


A- 1 


A 


A- 1 


A- 1 


A 


A- 1 


A 



(ni,n 2 ) 


Loops 




(0,0) 


1 


1 


(-1,0) 


U>4 + W^Ws + W4W7W8 + W3W4W7WS 
+Wi W5 Wq + W± 1 W§ 1 + W± + 1 


A 


(-2,0) 


Wi W5 W4 + W4W8 + W^W^Ws + Wi W$ W<±W8 
+W± W5 Wq W4W8 + W4W7WS + Wi W4W7W8 + Wi W§ W4W7W8 
+W3W4W7W8 + Wi W3W4W7W8 + W^ W§ W3W4W7W8 + W^Wq Wq 
+Wi 1 W§ 1 + W3WIW7W8 + WiW^W^WQ 1 + W^wlw^wl 


A 2 


(-3,0) 


W± W§ W4W8 + Wi 1 W^ 1 Wq 1 W 4 W8 + W^W^W^WfWs 

+Wi 1 W^ 1 W3WiW7W8 + W^Wq WZW4W7W& + W 3 w\w7w1 

+Wi 1 W^w\w7w'l + Wi 1 W$ 1 W 3 wlw7wl 


A 


(-4,0) 


w^ 1 w^ 1 w 3 wlw 7 wl 


1 


(-2,1) 


Wi 1 WiW7W8 


A 2 


(-2,-1) 


W2W3WIW7W8 


A 2 



(5.9) 
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Higgsing 2: Likewise, we tabulate the result for Higgsing 2 and obtain: 



Oi,n 2 ) 



Loops 



(0,0) 



1 



(-1,0) 



W4 + W4W8 + W4W7WS + W 3 W4W 7 Wg 



+w 1 1 w 5 1 w 6 1 + w 1 1 w 5 1 + w 1 1 



1 



W^ 1 W^ l W4 + WjWs + W^WiWs + W^W^W^Wg 
+W : [ 1 W^ 1 Wq 1 W4Ws + W4W7WS + W^W^W^Ws + W± W5 W4W7W8 
+ W3W4W7WS + W^WsWiWYWs + Wi 1 W$ W3W4W7W8 + Wi W§ Wq 



-2,0) 



-W 1 W 5 1 + W3WIW7W8 + W4W l W 5 W % 1 



w 3 wlw 7 wj 



-3,0) 



W^W^W^Ws + Wi W§ Wq W4WS + Wi W§ W4W7W8 

+Wi 1 W^ 1 W3W4W7W8 + W^W^WzwIwjWg + ti^f H>7wf 



-4,0) 



w 1 w 5 1 w 3 w'lw 7 wl 



(-2,1) 



W-y 1 W4W7Wg 



A" 



-2,-r 



w 2 w 3 wiw 7 w 8 



A" 



(5.10) 



Figure [X3| shows the amcebas corresponding to the scalings in ( |5.9| ) and ( p,10[ ) 
confirming they produced the desired splitting. 



y 




y 




-6-4-2 2 4 6 



-2 



Higgsing 1 Higgsing 2 

Figure 13: Amoeba plots for Higgsings 1 and 2 of Y '° at A = e 3 . 
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6. Combining Multiple Components 

In the previous sections, we have explained how to split integrable systems. Reversing 
the logic, we also understand how to glue them. We have identified a continuous 
parameter A that controls the distance between components of the spectral curve. This 
parameter manifests itself in the associated quivers as non-zero vevs and suppresses 
certain contributions in the integrable system. We can now proceed towards our goal 
of understanding the continuous limit of these systems. 

6.1 Combinatorics of a Large Number of Components 

Let us focus on the case in which we combine an infinite number of identical components 
E along a single direction, effectively generating a new continuous dimension. The 
amoeba projection suggests a natural way to approach the continuum: we consider all 
components equally separated in the amoeba and then send the number of components 
N contained in a finite interval of length L to infinity. More concretely, defining A = e a , 
we consider the limit 

N -> oo, a -)- 0, L = Na fixed. (6.1) 

This will give us a (1 + l)-dimensional integrable system from the continuous limit of 
an infinite number of (0 + l)-dimensional ones. In principle, it seems possible to do the 
same in both the x and y directions, generating a (2 + l)-dimensional integrable field 
theory in the process. 

Before studying this limit, let us investigate how the number of contributions to 
Hamiltonians behaves for large N. Following the dictionary in Section 0, this number 
corresponds to the multiplicity of perfect matchings associated to internal points in the 
toric diagram. 6 

Gluing N copies of an integrable system corresponds to considering a certain Z^ 
orbifold of the basic theory. In dimer model language, this corresponds to enlarging 
the unit cell by a factor N. While the number of points in the toric diagram grows 
linearly with N, their multiplicity grows much faster. Let us illustrate this growth in 
some explicit examples. 



6.1.1 Y 



JV.O 



The cone over Y ' is the Zjy orbifold of the conifold with toric diagram given by 



Figure WA . 



6 Casimirs are given by ratios of external points in the toric diagram. Hence, the same ideas apply 
independently to their numerator and denominator. 
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N/2+1 



Figure 14: Toric diagram for Y '° for even N. The red circle indicates the reference 
perfect matching and the green dots correspond to cycles with windings (—N/2 — 1, 1) and 
(—N/2 — 1,-1), which are fixed by the Casimirs. 



The integrable system for this geometry was determined in |15fl using the prescrip 



tion in [[13|], where it was identified with the iV-site relativistic periodic Toda chain. 
Our goal in this section is to investigate its behavior for large N. 



hH 



a ii o — ••• — ii i — I — 



d, 

i=l,...,N 



<l C^— i i Q — ••• — !► 



F3 



Ci 

even i 



Ci-l 
eveni 



Figure 15: A convenient set of cycles for Y N, ° with even N. The cycles of type c only exist 
for even i. 



For simplifity, let us focus on the case of even N. A similar analysis is possile 
for odd N. It turns out that the resulting integrable system is considerably simplified 
when considering the basis of cycles given in Figure [H| instead of using the standard 
w and z- variables. 7 In terms of this basis, the Hamiltonians become 



7 Figure [15] shows 2N cycles. The two additional cycles that are necessary to form a basis, are fixed 
by the Casimirs and hence not important in our discussion. They correspond to the green dots in 
Figure |TJ. 
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n factors 
The problem of finding the Hamiltonians is thus reduced to the combinatorics of non- 
intersecting paths, which can be used to immediately determine the multiplicity of 
internal points in the toric diagram. A closed expression for this multiplicity was 
derived in |26| by using a Potts model like description for the dimers, and via a recursion 
relation that was obtained from a map to a 1-dimensional monomer-dimer system. The 
final result for the multiplicity of the n th internal point is 



2-^ N -i 



j=0 



N-i 



n 



(6.3) 



which applies for both even and odd N. 

It is interesting to visualize how these multiplicities are distributed over the toric 
diagram and how the distribution approaches some limit shape after appropriate nor- 
malization. This is shown in Figure |16|. 



Mult. 




n/N 



Figure 16: Normalized multiplicity of perfect matchings for the internal points of the Y ,0 
theory. We have also normalized the length of the toric diagram to 1. We show results for 
N = 2a, a = 1, ... ,20 (black to red). 



6.1.2 Multiple F 



Similarly, we can investigate the generalization of the model considered in Section |5T 
to iV copies of Fq. Figure [17| shows its toric diagram. 

In this case, it is also possible to find closed formulas for the multiplicities of all 
points in the toric diagram. They are 
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• • • 



Figure 17: Toric diagram for the N Fq model. 



Boundary points: 



N 



m 



Internal points: 2 



2N 
2n-l 



(6.4) 



where N is the number of Fo's that have been glued together and m G {0,iV} and 
n G {1,N} index points on the boundary and in the interior of the toric diagram, 
respectively. 



Mult. 




Figure 18: Normalized multiplicity of perfect matchings for the internal points of the multi- 
ple Fq theory. We have also normalized the length of the toric diagram to 1. We show results 
for N = la + 1, a = 0, . . . , 20 (black to red). 



These two examples illustrate a general behavior of large- iV models, an explosive 
growth in the number of perfect matchings associated to a given point in the toric 
diagram, which translated into a huge number of contributions for each conserved 
charge. It thus becomes clear that a continuous reformulation of cluster integrable 
systems is desirable in order to deal with their large- N limit. In the next section we 
take the first steps towards such a reformulation. 
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6.2 A Toy Model for the Continuous Limit 

The two examples considered in the previous subsection share some common charac- 
teristics. In both of them, the n th Hamiltonian corresponds to the sum over all possible 
positions on the brane tiling of n paths, subject to the constraint of not overlapping 
over edges. Furthermore, these paths are of a very specific type: they cross the tiling 
along the short direction of the unit cell and are almost straight. By this we mean 
that these paths are almost localized along the long direction of the unit cell. Fig- 
15| shows the explicit form of these paths for Y N, ° and (|6.2|) gives the corresponding 



ure 

Hamiltonians. It is natural to assume that this structure is generic when gluing N 
copies of a cluster integrable system with a genus-1 spectral curve. 8 In this section, 
we introduce a toy model with these properties, which we expect captures the main 
features of a continuous reformulation of cluster integrable systems. 

Let us consider a system in the x £ [0, L] interval and introduce a path z winding 
vertically at x — 0. 



a 












t 


i 
z' 


W\ 


w 2 










w k 





x=0 x=L 

Figure 19: A toy model for the continuous limit of cluster integrable systems. The path z 
can be shifted by multiplying by all the Wj's contained in the strip between its initial and 
final positions. 



We can shift z horizontally by multiplying it by all the Wi variables contained in 
the strip between x = and its final position. 



z — >> 



J] W lZ . (6.5) 

iestrip 

Here Wj stands in general for a product of w^s contained in a slice of the tiling. 
The explicit form of Wj is controlled by the details of the specific brane tiling under 
consideration. Combining n paths and summing over their positions k iy i = 1, . . . ,n, 
we obtain 



8 Attaching N copies of a genus- N spectral curve results in gN Hamiltonians. 
resulting theories to obey a similar structure. 



We expect the 
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In ki \ 

yi=l ki<k %+1 j=0 J 

which we have suggestively called H n , since it has the expected structure for Hamil- 
tonian operators. 9 We can readily generalize this expression in the continuous limit, 
obtaining 

Hn=(f[[ L dx^ d y^v)\ A (6 . 7) 

where xo = 0. While we have derived QS.7| ) under rather basic assumptions that try 
to capture the most basic features observed in explicit examples, we expect it displays 
the main aspects of the actual continuous limit of cluster integrable systems. We leave 
a detailed investigation of this limit in explicit models for future work. 



7. Conclusions 

We have taken the initial steps in extending the correspondence between dimer mod- 
els and (0+l)-dimensional cluster integrable systems to continuous (1+1) and (2+1)- 
dimensional integrable theories. In order to understand the transition between discrete 
and continuous theories, it is necessary to have certain notion of distance between el- 
ementary constituents, or "lattice spacing", such that the continuous theory emerges 
when it is sent to zero. We identified such a continuous parameter controlling the 
distance between daughters from the perspectives of both spectral curves and the res- 
olution of Calabi-Yau singularities, equivalently the Higgsing of quivers. Furthermore, 
we introduced two procedures for determining the integrable system dependence on 
this parameter, whose effect is to suppress certain contributions to conserved charges, 
making them vanish in the infinite separation limit. 

We then explored the integrable systems that are constructed by combining a large 
number of components, equivalently by gluing a large number of toric diagrams. More 
concretely, we studied, in explicit examples, the behavior of the number of contributions 
to individual Hamiltonians as the number of components grows. These contributions 
are in one-to-one correspondence with perfect matchings of the underlying dimer model. 



9 Strictly speaking, Hamiltonians might also contain an n-independent power of the cycle orthogonal 
to z, as dictated by the position of the corresponding internal point in the toric diagram. This fact 
can be trivially incorporated in our expressions so, for simplicity, we omit it from our discussion. 
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For this reason, their number diverges much more rapidly than the number of compo- 
nents, begging for an alternative continuous formulation of cluster integrable systems. 
We also observed that each Hamiltonian is given by the contributions of a number of 
simple paths summed over all their possible positions on the brane tiling. We used 
these insights to develop a toy model that we expect reproduces the basic features of 
the continuous limit of cluster integrable systems. 

Interestingly, our investigation of the continuous merging of integrable system has 
also resulted in a novel understanding of (un)Higgsing in quiver theories and the asso- 
ciated desingularization of the corresponding Calabi-Yau spaces. Thus, we have added 
a new angle of attack to the classical subject of D-brane resolution of singularities. We 
have realized that when one refines the coefficients of the spectral curve into polyno- 
mials in loop variables, the Higgsing/resolution simply corresponds to establishing a 
consistent scale A dictating which monomials should survive or suppressed. Therefore, 
we have effectively generated a new algorithm, outlined in Section [4.2| , for systemati- 
cally studying all partial resolutions for a given toric diagram. It would be worthwhile 
to exploit this procedure for classifying all consistent daughter theories for a given 
parent. 

What we have touched upon is, of course, only the beginning of a program. The 
natural question that arises now is how to extend our continuous toy model to theories 
that are actually constructible from dimers. For example, attempting to recover simple 
integrable field theories such as Toda theories would be an obvious next step. Given the 
simplifications afforded by dimer models, we expect it should be possible to construct 
increasingly more elaborate integrable field theories. 

We envision many applications of dimer models to continuous theories, such as 
the study of integrablity preserving lower dimensional impurities or interfaces between 
different integrable field theories. Indeed, in light of the correspondence with dimer 
models, the often difficult condition of integrability simply amounts to checking whether 
the field theory results from the infinite limit of consistent toric diagrams, i.e. that they 
are given by convex lattice polygons. 

In addition to the continuous limit, there are several exciting directions worth 
studying for the dimer model/integrable system correspondence. For example, the cor- 
respondence naturally associates (0+l)-dimensional relativistic integrable systems to 
the (3+l)-dimensional superconformal field theories dual to the dimer models. Re- 
cently, similar integrable systems have emerged in the study of (3+l)-dimensional 



M = 2 superconformal theories in the contexts of superconformal indices [27j and 
the enumeration of vacua in the Omega background [^, |29| . It would be interesting to 
investigate whether the different ways in which integrable systems emerge are indeed 
related. 
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